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Abstract 



1 Introduction 



The radiative decay B — > X s 7 plays an important role in the present tests of the Standard 
Model (SM) and of its extensions [0. In particular, in the supersymmetric extensions of 
the SM, the best bounds on several new parameters come from the data on B — > 0] . 

The short distance QCD effects are very important for this decay. They are known to 
enhance the branching ratio Br(5 — > X^) by roughly a factor of three, as first pointed 
out in Since these first analyses, a lot of progress has been made in calculating 

these important QCD effects in the renormalization group improved perturbation theory, 
beginning with the work in || Let us briefly summarize this progress. 

A peculiar feature of the renormalization group analysis in B —>■ Xgj is that the 
mixing under infinite renormalization between the four-fermion operators {Qi, ■ ■ ■ ,Qe) 
and the magnetic penguin operators (Q7 7 , Qsg), which govern this decay, vanishes at the 
one-loop level. Consequently, in order to calculate the coefficients CV 7 (/ij,) and Cscif^b) 
at /lb = 0{m h ) in the leading logarithmic approximation, two-loop calculations of 0(eg1) 
and O(g^) are necessary. Such calculations were completed in |7], [| and confirmed in 
[H [lO], |llj. Earlier analyses contained either additional approximations or mistakes. 

It turns out that the leading order expression for the branching ratio Br(5 — > X^) 



suffers from sizable renormalization scale uncertainties |12|, [L3| implying that a complete 
NLO analysis including also dominant higher order electroweak effects to this decay is 
mandatory. By 1998, the main ingredients of such an analysis had been calculated. It 
was a joint effort of many groups: 



The 0(a s ) corrections to Cy^nw) and Cgc(nw) were first calculated in [jTj] and 
subsequently confirmed by several groups fT5||-fli 



The two-loop mixing involving the four fermion operators Qi, ■ ■ ■ ,Qe and the two- 
loop mixing in the sector {Q^, Qsg) was calculated in |19||- [^| and fl25f , respectively. 
The very difficult three-loop mixing between the set (Qi, . . . , Q 6 ) and the operators 
(Q77) Qsg) was calculated in [[U] . 



Leading order matrix elements (sj gluon|Qi|&) were calculated in ^7|, and the 
challenging two-loop calculation of (s r y\Q2\b) was presented in EE . 



Higher order electroweak corrections were incorporated with increasing level of so- 



phistication in p9|-p3 



In addition, non-perturbative corrections were calculated in [S3]-[5S|. The most recent 



analysis of B — > X s ^ incorporating all these calculations can be found in [^0 . 
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Now, among the perturbative ingredients listed above, three have been calculated only 
by one group. These are 

• The two-loop mixing in the sector {Qt 1: Qsg) HH- 

• The three-loop mixing between the set (Qi, . . . ,Qq) and the operators (<57 7 ,<58g) 

• The two- loop matrix element (sj\Q2\b) |p8| . 
Moreover, 

• The two-loop matrix elements of the QCD penguin operators (s7|Qj|6) with 
i — 3, ... ,6 have not yet been calculated. 

It should be emphasized that all these four ingredients enter not only the analysis of 
B — > Xg'-f in the SM but are also necessary ingredients of any analysis of this decay in 
the extensions of this model. It is therefore desirable to check the first three calculations 
and to perform the last one. 

In the present paper, we will make the first step in this direction by calculating the two- 
loop matrix element (s r y\Q2\b) using the method of asymptotic expansions. This matrix 
element turned out |^8|, to be the most important ingredient of the NLO-analysis 
for Bt(B — > Xg'-f) enhancing this branching ratio by roughly 20%. In ||28| , the matrix 
element (S7IQ2I&) wa s found by applying the Mellin-Barnes representation to certain 
internal propagators, and the result was presented as an expansion in z = m 2 c /m1 up to 
and including terms 0(z 3 ). In order to be sure about the convergence of this expansion, 
we will include also the terms 0(z A ), 0(z 5 ) and 0(z 6 ). 



2 Two-Loop Contribution to (s r y\Q2\b). 
2.1 Preface 

The effective Hamiltonian for the process b — > 57 is given by 

AC 8 

H eS (b -> s 7 ) = ^ Vth V t lY,Ci^)Qi- (2-1) 

Here Gp is the Fermi constant, Vij are the CKM matrix elements and Cj are the Wilson 
coefficients of the operators Qi evaluated at /i = 0{m b ). We have dropped the negligible 
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contributions proportional to V u bV* s . The full list of operators in ( |2.1[ ) can be found in 
[ 28| . In the present paper we need only the expressions for two of them. These are 

^2 = (c^P L b)(s^P L b), (2.2) 
Qi = -^— 2 s a o^(m b (ii)P R + m s (ii)P L )b a F^, (2.3) 



whereP L = (1-J 5 )/2,P R = (l+7 5 )/2 and = f[7 M ,7„]. As in |28| , we will set m s = 



In this section, we present the details of the calculation of the matrix element (s , y\Q2\b) 
in the NDR scheme. This matrix element vanishes at the one-loop level. Therefore, in 
order to obtain a non-vanishing contribution, one has to calculate two-loop diagrams. 
They are shown in Figs. |], f|, ||, and || where the wavy and dashed lines represent the 
photon and gluon, respectively. Following f2~8fl , we have divided the contributing diagrams 
into four sets: 

• In sets 1 and 2, the photon is emitted from an internal s-quark (setl) or an internal 
b-quark (set 2), which also exchanges a gluon with the charm quark. 

• Sets 3 and 4 are obtained from sets 1 and 2, respectively, by emitting this time the 
photon from the charm quark propagator. 

• Thanks to QED gauge invariance, it is not necessary to consider diagrams with a 
real photon emission from external quark lines. 

It is convenient to write the regularized contributions from each set of the diagrams 
as follows: 

M l = loirf— ) +M 1 )-C F Q d (si\Q 7 \b)treei (2.4) 
I 306 \K[1 J I 7T 

M-2 = {-£r(—) +^1-^,(57^716}^, (2.5) 

3oe J 7T 

M :i = {~(^\ +M 3 \^C F Q u ( S1 \Q 7 \b) tTee , (2.6) 

8e \Kfl J I 77 

7/, = {-^( T ^) + MA^C F Q u ( S1 \Q 7 \b) tiee . (2.7) 
4e \KfM J 7T 

We work in D = 4 — 2e dimensions, with k 2 = 47re -7B , Q u = |, = — |, C F = | and 
je is Euler's constant. The tree level matrix element of the operator Qj is given by 

g 

(s7|Q 7 |6)tree = m 6 (// 6 ) —u(p')jt<£P R u(p) . (2.8) 
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In order to make the two-loop matrix element finite, counterterms have to be added. 
These counterterms can easily be calculated by using the operator renormalization con- 
stants that were needed in the context of the calculation of the leading order anomalous 
dimension matrix. The complete counterterm is found to be 28| 



A/7 xount \Qd(m b \ 3Q U Q d [ n 4t a s 

M 2 =\j[—) + — -^j--C F (si\Q 7 \b) tme . (2.9) 

The above expression includes two-loop counterterms as well as contributions from one- 
loop diagrams with one-loop counterterm insertions. 

Adding the contributions from (|2.4|) - (|2.7|) and (|2~^ ), we find the two-loop matrix 
element (sj\Q2\b) in the NDR-scheme: 

( S7 |Q 2 | fe ) = ( S7 |Q 7 | 6)trec ^ ^ hi V± + r^j , (2.10) 

with 

r 2 = - — (Mi + M 2 ) + —(M 3 + M 4 ). (2.11) 
y » 

The rest of this section is devoted to a detailed presentation of the calculation of the 
contributions Mj using the method of asymptotic expansions. All the momentum integrals 
are performed in the Euclidean space. 



2.2 Diagram M x 

Of the four sets of diagrams which we have to consider, Mi and M 2 shown in Fig. [I] and 
H are relatively simple. In M 1; the photon is emitted from a massless s-quark propagator, 
and the two-loop integral factorizes into the c-quark loop and the one-loop vertex integral. 
There are four momentum regions we have to consider. Denoting the momenta in the 
c-quark and vertex loops by p and k, we have 

1. k ~ p ~ nib: the "hard-hard" region; 

2. k ~ mb andp ~ m c : the "hard-soft" region (which comes in two variants, depending 
on the routing of k through the c-quark loop); 

3. a "collinear-soft" region, in which p ~ m c and k ■ 7 ~ ml but k ■ q ~ m\. 
Below, we describe the treatment of those regions in some detail. 
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Figure 1: Diagram Mi. Momentum assignments: Q, q, and 7 denote respectively four- 
momenta of the b and s quarks, and of the photon; p and k are loop momenta. Double 
lines denote a massive particle (the b quark) and single lines are light particle propagators. 
The gluon is denoted by a dashed line. 



2.2.1 Hard contribution: Taylor expansion in m c 

Since all scalar products are now large compared with m 2 , we can expand the c-quark 
propagators in m c . This leads to a massless loop integration, which simply modifies 
the power of the momentum (k) in the gluon propagator (see Fig. |2|). After this first 
integration, the integrand has the following form: 

d D k (scalar products involving k) 



( fc2 )a 1+ e(p + 2k-Q + Q 2 )^(k 2 + 2k ■ q)^' ^ 2 ' 12 ^ 

We combine propagators 1 and 3 using the Feynman parameter x, and then combine the 
result (-(1 — y)) with propagator 2 (-y). We shift the variable k K — x(l — y)q — yQ, 
and the denominator becomes a power of k 2 + y(l — y)(x — 1 — iO). Since x — 1 < 0, 
this diagram has an imaginary part; its sign is determined by assigning a small negative 
imaginary mass to the massless lines. After the momentum shift, we have additional 
factors x m and y n in the numerator, and denominator has the power d + e. The result is 

2 d r(d-2 + 2e) B(m + a 3 ,a 1 + 2-d-e) 



|exp(— in)] e (— 1 , 

L J T{d + e) B(a 3 ,a 1 + e)B(a 2 ,a 1 + a 3 + e) 

xB(n + a 2 + 2 - d - 2e, a x + a 3 + 2 - d - e), (2.13) 

where B and V are the standard Euler functions Beta and V. 



2.2.2 Hard-soft contributions from k 



m b , p ~ m c 



This region can be handled as in the "large momentum expansion" [[IT], 42, 43]. In Fig. || 
we see one of the two configurations to be computed. We first integrate over p, and then 
evaluate the one-loop vertex. An analogous procedure is executed for the other (left) 
c-quark line. 
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Figure 2: One-loop diagram (the k propagator can have fractional power) corresponding 
to eq. (|212p . 



\k 



q+k 





Q " Q+ k \ q 

y 



Figure 3: Subgraph corresponding to one of the two hard-soft contributions. 



2.2.3 Collinear region 

Naively, one might think that we should consider the region k ~ p ~ m c , in which 
(Q + k) 2 ~ Q 2 , (q + k) 2 ~ 2k • q. However, this does not lead to a consistent power 
counting, since we have to shift the momentum k by an amount proportional to q, which 
is ~ m&. 



Instead, we have to consider the collinear configurations [44|. We define 

k-q = Q k_, k-'y = Q k + , k 2 = k + k_ - k 2 ± , Q 2 = Q 2 Q) k ■ Q = Q (k + + k_). 

(2.14) 

2 

The relevant contribution is k + ~ Qq, k- ~ Now (Q + k) 2 ~ Q 2 + 2k ■ 7, while 
(fc + q) 2 = k 2 + 2fc ■ q is already homogeneously of order and cannot be expanded. 
First, we integrate over p; for this we combine the two c propagators (we assume here, 
for generality, that their powers are ai and 02) using x, shift p ^ P — xk, and obtain (m 
and d are exponents arising from the presence of p in the numerator) 

polynomial^) T(d - 2 + e) afn+ai+i-d-e^ _ x )«2+i-d-e 



[(p + A;) 2 + mf 1 [p 2 + m^f 2 S(a 1 ,a 2 )r(d) (k 2 + m 2 x ) d - 2+t 



> m 2 

m 2 x = c 2.15 

x(l - x) 
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For the k integration, the integrand has the form 

d D k polynomial^) 



( fc 2)a 3 (fc2 + m 2y-2+e( k 2 + ^ . g)«4 (Q2 + 2? . ' 



(2.16) 



Now, we combine the first two propagators multiplying them by 1 — z and z, respectively. 
We multiply the result by y and combine it with the "a 4 " propagator multiplied by 1 — y. 
Finally, we use the parameter u to include the last propagator (i = a 3 + a 4 + d — 2 + e) : 

1 



[k 2 + yzm 2 + 2(1 - y)k ■ q\ l (Q 2 + 2k ■ 7)^ 

,05-1 



" B(i,a 5 ) Jo dU [k 2 + yzm 2 + 2(1 - y)k ■ q + 2mA; • 7 + wQ 2 ]^ ' 1 2 ' J ' ' 

Now, we shift k — > K — (1 — y)g — M7 and, using 5 • 7 = Q 2 /2, find that the denominator 
simplifies to become [i^ 2 +?/(zm^ + -uQ 2 )]. Because of the presence of in the numerator, 
we generate additional powers of Feynman parameters. After integrating over K, we 
integrate over u, 



f 

Jo 



00 u n 1 B(n,m — n) 

{zm 2 + uQ 2 ) m = {Q 2 ) n (zm 2 ) m - n ' ( ' ' 



We see that the dependence on Q separates (so the integral really depends only on a single 
scale). Also, only m c appears in fractional power, so that there is no imaginary part in 
this integral, even though Q 2 < (after the Wick rotation). 



2.3 M 2 




Figure 4: Diagram M 2 . 



In this diagram, we have again four momentum regions, of which the first three are 
similar to those discussed for Mi. The fourth one is different: it is a "soft-soft" contribu- 
tion, with k ~ p ~ m c . 
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For the hard-hard contribution, we first integrate over p, just like in M\. After that, 
the integrand takes on the following form: 

d D k (scalar products involving k) 

(2.19) 



( k 2y 1+ e( k 2 + 2k ■ Q)^(k 2 + 2k-q- Q 2 ) a * 

We combine propagators 2 and 3 using a Feynman parameter x, and then combine the 
result (•?/) with the propagator 1 (-(1— y)). We shift the variable k — > K — y((l — x)Q+xq), 
and the denominator becomes a power of k 2 +y[y+x(l — y)]. After shifting the momentum 
and simplifying the numerator, additional factors x m and y n appear in the numerator, 
and the denominator has the power d + 2e. After the k integration, we find that the x 
and y integrations do not separate, and we first integrate over x using eq. (|2.30| ), to be 
discussed below. The remaining y integration is a Beta function. 

The hard-soft contributions are similar to those in M\ and lead to rather trivial prod- 
ucts of one-loop integrals. 

On the other hand, the soft-soft contribution is less standard. We have k ~ m c and 
the two 6-quark propagators can be expanded. In the lowest order we have 

1 1 



k 2 + 2k-Q 2k -Q' 
1 1 



k 2 + 2k-q-Q 2 Q 2 ' 



[2.20) 



Since the second propagator becomes independent of k, the whole two-loop diagram be- 
comes equivalent to a simpler, two-point function. The resulting integrals are similar, but 
not the same, as those considered in the eikonal expansion study P5 |. Details of their 



evaluation have been recently described in a completely different context of bound-state 



calculations | 46fl . In that work, exactly such type of integrals appeared when energies 
of bound-states consisting of two particles with very different masses were expressed as 
expansions in the ratio of those masses. 

2.4 M 3 

This diagram, shown in Fig. |^, is the first of the two non-trivial two-loop vertex diagrams 
we have to consider. There are now five momentum regions to be considered: 

1. Hard-hard, k ~ p ~ mj, similar to those of M 12 . 

2. Hard-soft, k ~ mj, p ~ m c , which now enters only once, when hard momentum 
flows through the same c-quark line from which the photon is emitted. 
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k+q q q k+q q 

(a) (b) 
Figure 5: Diagrams contributing to M3, with photon emission from two c-quark lines. 

3. Collinear-collinear, with k and p having their only large components [O aligned 
with the s-quark momentum q (and with k 2 ~ p 2 ~ m 2 ). 

4. Collinear-collinear, but with the alignment with the photon momentum 7. 

5. Hard-collinear, where k ~ mj and p is aligned parallel with 7. 

The first two regions are treated in an analogous manner as was described for M-yp- 
In the hard-hard contribution, after expansion in m c , the integrand is symmetrical (apart 
from different powers of propagators) under replacement p <-> k, 7 «-> g. Hence, we can 
first integrate over k, and then over p. Both integrations are very similar to those in Mi 2. 



In the second region, we use the "large momentum expansion" mentioned in Section |2.2.2 . 

The collinear-collinear regions require more attention [|47|| . After expansion of the in- 
tegrands in the available small quantities, it turns out that the integrals over Feynman 
parameters have singularities (like J dx/x) which are not regularized by our dimensional 
regulator. It is necessary to introduce additional, analytical regularization on the heavy 



quark line. Similar phenomena have been observed before (see e.g. [57| and references 



therein). Nevertheless, all integrals over the Feynman parameters can be evaluated an- 
alytically without particular difficulties. Dependence on the analytical regulator cancels 
in the sum of the two doubly-collinear contributions. 

2.5 M 4 

The most complicated diagram is M4, depicted in Fig. [f| There are six momentum regions: 
1. Hard-hard. 



2. Hard-soft. 



3. Collinear-collinear with alignment along 7. There is no contribution with alignment 
along q here, because of the different structure of the internal quark propagator, 
which now contains the large 6-quark mass. 



4. Hard-collinear, as in M 3 . 

5. Ultrasoft-collinear [I7|, with k ~ ml/m^ and p aligned with 7. 

6. Soft-soft. 





Q Q+k s Q Q+k s 

(a) (b) 
Figure 6: Diagram M4 

We give a more detailed description of the hard-hard contribution, because the inte- 
grals resulting here are rather complicated. First, we consider the scalar integral with all 
propagators in the first power, 

d D k r d D p 



Js J k 2 (k 2 + 2Q-k)J p 2 (p-k) 2 (p + 1 ) 2 ' {2 ' 21 ^ 

We first evaluate the massless (c quark) loop (integrals over all Feynman parameters run 
from to 1): 



d D p 



p 2 (p — k) 2 (p + 7) 5 



d D p 



p 2 (p 2 — 2p ■ k + k 2 )(p 2 + 2p ■ 7) 



/ 2(l-y)dydx J — 



d u p 



[p 2 — 2yp ■ k + 2x(l — y)p ■ 7 + yk 



213 



r(l + eW 



2-e 



' l e (l — y) e dy dx 



(k 2 + 2xk ■ 7) 



l+e 



T(l 



)B(-eA-e)7i 



2-e 



dx 



(2.22) 



(k 2 + 2xk- 1 ) 1 +'' 

In the process of the integration, we made a shift p — > p + yk — x(l — y)j. Next, we 
integrate over k: 

d D k 



dx 



k 2 {k 2 + 2Q-k){k 2 + 2xk- 1 ) 1 +" 
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;i + e) j dx du (l-u) e J 



d D k 



k 2 [k 2 + 2uk-Q + 2(1 - u)xk ■ 7 ] 2+e 
(2 + e)(l + e) J dz z 1+e dx du (1 - uf J ■ <! ' ''' 



[k 2 + z 2 u 2 + z 2 u(l — u)x] 3+e 



r(l + 2e) 2 _ e f. . . z- 1 - 3 ^- 1 - 26 ^ - u 
V e / dz dx du 



r(l + e) J [u + (1 -u)x]i+2« 

1,1 7T 2 - £ /"dtttt-^Cl-^-^Cl-O 



6e 2 r(l + e) 
T(l + 2e) 
6e 2 r(l + e) 



-2-e 



7T 



[B(-2e,e)-B(-4e,e)}. (2.23) 



We made a shift k — > — zuQ — ^(1 — m)x7 and used Q 2 = 2Q ■ 7 = — 1. Multiplying the 
result (12T23D with the coefficient T(l + e)S(-e, 1 - e) from (ggp we find 



J 3 = V ,7 2 1 B(-e, 1 - e) [B(-2e, e) - B(-4e, e)] 



r(l_+2e) 

6e 2 

= r(l + e ) 2 f^- + ^- + -C 3 + ^V^(e)- (2.24) 
v ; \24e 4 18e 2 6e S 360 J w V 7 

Now we consider the general case, in which momenta p and k can be present in the 
numerator, 



d k d p (scalar products involving k and p) 



(2.25) 



{k 2 ) ai (k 2 + 2k ■ Qf 2 (p - k) 2az (p 2 ) ai {p 2 + 2p ■ 7 ) as ' 

We first combine the propagators 5 (-x) and 4 (-(1 — x)), and then the result (-(1 — y)) 
with 3 (-y). We shift the integration momentum p = K — x(l — y)^y + yk and average over 
K. This again results in extra powers y n and x m , while canceling powers of K 2 changes 
the power of the denominator. After the integration over K, the denominator becomes 
[y(l — y)(k 2 + 2xk ■ 7)] d_2+<E . We see that the dependence on y factorizes and we can 
integrate over this variable. As a result, we find an expression of the form 

T(a 3 + a 4 + a 5 )T(d - 2 + e)B(n + a 3 - d + 2 - e, a 4 + a 5 - d + 2 - e) 

r(a 3 )r(a 4 )r(a 5 )r(rf) ( ' j 

d« a5+m - 1 (l-i) a4 - 1 (2.27) 



/" d D k (scalar products involving k) 

x / w 7 . 2.28 

J {k 2 ) ai {k 2 + 2k-Q) a \k 2 + 2xk- 1 ) d ~ 2+e V ; 

We now repeat a similar procedure with the variable p. We combine the propagator 2 (-it) 
with (/c 2 + 2x/c -7) (-(1 — it)), and then the result (-z) with the propagator 1 (-(1 — z)). We 
change the momentum variable k = K — zuQ — zx(l — u) 7 , average over K and simplify 
powers of K 2 (as a result, the power of the denominator changes from ai + a 2 + d — 2 
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to some d\, and we also get extra powers z e and u w ). After integrating over K we get 
(including the factors ( ggg , gg7| )) 

T(a 3 + a 4 + a 5 )r(d — 2 + e)B{n + a 3 — d + 2 — e, a 4 + a 5 — d + 2 — e) 

r(a 3 )r(a 4 )r(a 5 )r(d) 

r(di - 2 + 2e)r(ai + a 2 + d - 2 + e)B(a u a 2 + d + 2 + e - 2d x - 3e) 



r(di + e)r(oi)r(a 2 )r(d - 2 + e 

dx d« x <l5+m_1 '' 1 — ^, N \a4-l„,«'+n2+l-rfi-2e / ' 1 _ n ,\d-3+e 



1 - x )«4-l^+a 2 +l-d 1 -26( 1 _ u y~-A+e ^ + ^ _ u )] 



. (2.29) 



First we integrate over x, using 



dx x n [it + x(l - it)]' 



u 



H+n+l n 



U 



in+1 



E( 

i=o 



n\ u 



j-X-n-fi 



J I fi + n + 1 - j 



[2.30) 



This integration is possible because the powers of x and (1 — x) in ( [2.271) are non-negative 
integer. The it integration gives a simple Beta function. 

The techniques described above are sufficient to compute all the remaining contribu- 
tions. Again, the dimensional regularization alone is insufficient to evaluate the doubly- 
collinear contribution. We introduce an analytical regulator on the heavy quark mass. 
The resulting singularities vanish when we add the soft-soft contribution. 



2.6 Results 

The final results for are listed below, with z = m 2 /m 2 and L = lnz. 
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L-h> 
4 



1 

12 



L 3 ] z 2 



+ 7T Z 

A 4 
'58 5 

+ I — - — 7T 



+ 



+ 



27 1* 
Z6547 _ 35 2 

V1728 ~ 72^ 
185837 21 



41 
72" 



2 --L 



L 2 )z i 



283 



144 



35 
24 



L 2 )z 5 



24000 20^ 



4227^ 63^2 
~800~ ~~ 20 



z . 



(2.34) 



For the imaginary part (and the leading power of L in the real part) it is possible to 
guess the form of the higher order terms. For example, for M 1; beginning with z A we have 



Higher orders in Im Mi = 8^ 3 



(2n + l)! 



{n (n + 1)! (n + 3)!' 



(2.35) 



If m c = rrib/2, we expect the imaginary part to vanish. Indeed, this corresponds to 
z — 1/4, at which point fl2.35|) gives ^jg — ^ In 2, which exactly cancels the contribution of 
the first 4 terms in the imaginary part of Mi given in ( |2.31j) . The results for the leading 
terms and the powers z, z 3//2 , z 2 and z 3 in ( j2.31| )-( |2.34| ) agree with the corresponding 
terms in (2.25)-(2.28) of Greub, Hurth and Wyler |28]. The contributions z 4 , z 5 and z 6 
are new. 

Inserting (g3TJ) (g34[) into (PUD we find 



T2 = Re T2 + 21H1 T2 



with 
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Rer 2 



Im r 2 



1666 32 9 , /9 

h — 7rV /2 
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+ 
+ 
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+ 

+ 



243 
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27 

— TT L -71 

3 27 



32 > 32 r 
3 9 27 9 27 



-L 2 
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4 56 

7T* + 

9 81 



728 



56 



-L L" 
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111748 


176 2 


3308^ 
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6075 


81 


405 


816731 


380 2 

7T - 


13234 


23814 


81 


567 



176 

27 
380 



/44551813 1624 



-TT 



421121 
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L — 

V 546750 135 6075 

/16 r 16 r2 80 16 ? 
L + —L + — - —tt j 



81 



9 



9 



64 448 

L H 

27 81 



16 



z° + — z 1 + 
81 



9 27 
100 



1624 
^5~ J 

z + 



— T 3 



z 



16 2 16 , 

L 7T 

9 27 



-z° + 



5516 
2025' 



;2.36) 



[2.37) 



confirming the results (2.37) and (2.38) of PB[ , and generalizing them to include terms 
0(z 4 ), 0(z 5 ) and to 0(z 6 ). The contributions of different terms z n for m c /m;, = 0.29 and 
m c /rrib = 0.22 are presented in table [IJ The term z 3 ^ 2 has been added to the term z 2 in 
this table. 





Rer 2 


Im 


r 2 


n 


2* = 0.29 


^ = 0.22 

m h 


^ = 0.29 







-6.8559671 


-6.8559671 


-3.1028076 


-3.1028076 


1 


2.2271232 


1.6504869 


2.5193530 


2.1225820 


2 


0.5775520 


0.2307186 


0.1121646 


0.0769360 


3 


-0.0402440 


-0.0103793 


0.0213018 


0.0045268 


4 


-0.0011397 


-0.0002090 


0.0000310 


0.0000034 


5 


-0.0001703 


-0.0000187 


0.0000163 


0.0000010 


6 


-0.0000307 


-0.0000020 


0.0000030 


0.0000001 



Table 1: The numerical value of Coefficient -z n for the real and imaginary part of r 2 are 
given up to z 6 using two different values of m c /mb. The term z 3 ^ 2 in Rer 2 was included 
in the term proportional to z 2 . 



We observe that the terms 0(z n ) with n > 4 are negligible. The final result for r 2 is 
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given by: 
Rer 2 = 



-4.093 
-4.985 



m c /m b 
m c /m b 



0.29 
0.22 



Im r 2 



-0.450 
-0.899 



m c /m b 
m c /m b 



0.29 
0.22 



(2.38) 

The strong dependence of r 2 on m c /m b has been pointed out in j5J. With decreasing 
z the enhancement of B — > X s j by QCD corrections becomes stronger. 



3 Conclusion 

In the present paper, we have calculated the important two- loop matrix element (s r y\Q2\b) 
contributing to the decay B — > X s ^ at the NLO level. Our result for (sj\Q 2 \b) agrees 
with the one presented in p8 and used in the literature by many authors. The additional 



terms in the expansion in z = m 2 c /m b , that is 0(z 4 ) and higher, turn out to amount to 
at most 0.05% and are negligible. As we have used a completely different method from 



the one used in [ 28] , the confirmation of the result of these authors is very gratifying. 
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